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Abstract 

We present a detailed study of a Brownian particle driven by Carnot-type 
refrigerating protocol operating between two thermal baths. Both the un¬ 
derdamped as well as the overdamped limits are investigated. The particle 
is in a harmonic potential with time-periodic strength that drives the system 
cyclically between the baths. Each cycle consists of two isothermal steps 
at different temperatures and two adiabatic steps connecting them. Besides 
working as a stochastic refrigerator, it is shown analytically that in the qua¬ 
sistatic regime the system can also act as stochastic heater, depending on the 
bath temperatures. Interestingly, in non-quasistatic regime, our system can 
even work as a stochastic heat engine for certain range of cycle time and bath 
temperatures. We show that the operation of this engine is not reliable. The 
fluctuations of stochastic efficiency/coefficient of performance (COP) domi¬ 
nate their mean values. Their distributions show power law tails, however the 
exponents are not universal. Our study reveals that microscopic machines 
are not the microscopic equivalent of the macroscopic machines that we come 
across in our daily life. We find that there is no one to one correspondence 
between the performance of our system under engine protocol and its reverse. 
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1. Introduction 


Thermodynamics of micro- and nano-scale systems exhibits distinctly dif¬ 
ferent features from that of large systems due to influence of large thermal 
fluctuations [Tj. Typical energy changes are of the order of thermal energy 
per degree of freedom and consequently thermodynamics has to be modified 
at micro-scale. These systems can be theoretically analyzed using stochas¬ 
tic thermodynamics. Exchange of energy between the particle and its sur¬ 
roundings becomes stochastic and yet one can clearly formulate the notion 
of work, heat and entropy production for a given microscopic trajectory of 
the particle [a El m El EU3 u E]. Recently obtained exact results (fluctua¬ 
tion theorems[TTj flOl 12, lEJ Q, [15]) put constraints on the distributions of 
the above mentioned stochastic quantities and are valid for systems driven 
far from equilibrium. These theorems transform thermodynamic inequalities 
into equalities. This area has become even more interesting with the develop¬ 
ment of experimental techniques. Using single-colloidal particle experiments 
several new key concepts have been verified. Information to energy conver¬ 
sion and validation of generalized Jarzynski equality [T6], Landauer erasure 
principle ira. universal features in the energetics of symmetry breaking [18] 
are to name a few. Micron sized heat engines have been experimentally 
realized by optically controlled motion of trapped colloidal particle HSH2Q]. 

There are several extensive studies on single bath nano-machines e.g., 
information machines (that can produce work using available information) 
|21, 22] and molecular motors / thermal ratchets [23, [24]. Molecular motors 
are omnipresent in cellular as well as tissue level of many living organisms. 
They are efficient enough to extract energy from a highly fluctuating environ¬ 
ment and to convert it into mechanical work for cellular and / or intra-cellular 
logistics [25] . 

The conversion of energy into mechanical work repeatedly along a cycle 
working within multiple thermal bath is almost every where in our day-to- 
day life, spanning huge range of length and time scales. For example, starting 
from high pressure steam locomotives, to a drinking bird [26] and even bio¬ 
chemical reaction pathways for cellular respiration mechanism [271128 J pro¬ 
ducing useful energy from nutrients - everywhere energy is transformed into 
mechanical work. Similarly, by reversing the cycle, we see that in various 
processes mechanical work is used to transfer heat from a cold source to a 
warm sink with an objective to cool down the cold source further (refriger¬ 
ator) or to heat up the warm sink warmer (heat pump). Though it is very 
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important to study the work-energy (or, vice versa) conversion in all relevant 
scales, due to lack of experimental techniques for micro or nano world, it is 
relatively well explored in macro scale. 

Heat engines and refrigerators at nano-scale is a subject of current study 
[2U1 m Ell E21 Em EH ESI ESI EZ] ■ Detailed theoretical treatment of Carnot- 
type micro heat engine, involving both quasistatic and non-quasistatic (i.e., 
finite cycle time) features, have also been documented [S3, 53 1 155 1 136; : 57l 138] . 
These features reveal the fundamental differences between micro and macro 
heat engines due to thermal fluctuations, reflected in the distributions of var¬ 
ious thermodynamic quantities (e.g., work, heat exchange, efficiency etc.). 
Unlike macro heat engines, it has been shown that the system can work as 
a heat engine if the ratio of hot and cold bath temperatures is larger than a 
critical value [38]. Moreover, in non-quasistatic regime, the system works as 
a heat engine for cycle times larger than a critical value. Both the thresholds 
depend on the system parameters. Fluctuations in thermodynamic quan¬ 
tities including efficiency of the system calculated over a large number of 
trajectories are significant not only in non-quasistatic regime but also in 
quasistatic regime, which is in clear contrast to the macro engines. Several 
trajectories violate typical expectations from the second law of thermody¬ 
namics [39] m\ ■ The non self-averaging nature of fluctuations in stochastic 
efficiency and other quantities requires detailed understanding of full proba¬ 
bility distributions as opposed to the average behavior mm- Large devi¬ 
ation properties of such distributions are recently under theoretical investi¬ 
gations Haim im ssi ii9ii5oi e]. Research on fluctuation relations for heat 
engines [52] [53, 54j are being pursued. It may also be noted that, for some 
of the heat engines studied so far, one may or may not recover Carnot result 
in the quasistatic regime. However, fluctuation theorems provide a bound on 
efficiency of an engine valid for any finite time cycle. Recently, novel theoret¬ 
ical approaches to capture the statistical properties of stochastic efficiency 
of micro heat engines and mesoscopic thermoelectric engines with broken 
time-reversal symmetry are being developed, particularly at long time limit 
[43, .H], claiming universal properties of the large deviation function related 
to the statistics of stochastic efficiency. To our knowledge, so far, no such 
study exists for micro-refrigerators. 

In this article we will focus on Carnot-type single-particle refrigerator 
and its stochastic features. The refrigeration protocol used here is similar 
to the micro heat engine protocol used in [38] but is running backward in 
time. We believe that our model system is experimentally realizable using the 
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technique already being used for micro heat engines. We find new insights 
into far from equilibrium features of the concerned system. For example, 
a major outcome of the present study is the variety of different modes of 
operation for such systems under the protocol. In Fig. [l]we describe all the 
modes of operations which are thermodynamically possible for a system that 
works cyclically between a hot and a cold heat bath. Other four possibilities 
for heat exchanges and work are ruled out due to violation of First and 
Second laws of thermodynamics. 



Figure 1: (Color online) Four thermodynamically possible machines working between a 
hot (with temperature Th) and a cold (with temperature 7]) thermal baths: A: Takes 
heat Qi from hot bath and converts it partially into work W and supplies the rest in the 
form of heat Q 2 to the cold bath (heat engine). B: Takes heat Q 1 from the hot bath and 
with the help of work W on it, supplies heat Q 2 to the cold bath. C: Converts work W 
on it, to heats Q\ and Q 2 that enter into the hot and cold baths respectively.B and C are 
called heater of type-II (or, heater-II) and heater of type-I (or, heater-I) respectively. D: 
Takes heat Q 2 from the cold bath and with the help of work W on it, supplies heat Q\ 
to the hot bath (refrigerator). When work is being done on (by) the system, it is positive 
(negative). When the system releases heat it is positive and negative otherwise. 


Our system, under the refrigeration protocol, exhibits all the modes de¬ 
pending on its cycle time and on the ratio of hot and cold bath temperatures, 
which can be controlled in experiments [19, 20j. Throughout the paper, we 
consider the heat going from the system to the bath to be positive and that 
going from the bath to the system to be negative. Similarly, the work done 
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on the system is taken as positive and that done by the system is considered 
negative. We have also studied fluctuation in stochastic efficiency / coef¬ 
ficient of performance (COP) of engine / refrigerator and have shown that 
their fluctuations dominate the mean values. Probability distributions have 
been calculated both in the underdamped and overdamped regimes. These 
distributions exhibit power law tails with non-universal exponents. Reliabil¬ 
ity of our system operating as a refrigerator / engine is also reported. 

We first describe the model and the protocol in the next section. In 
section 3 we discuss the essentials of stochastic thermodynamics. In section 
4 and 5, we explain the results for quasistatic and non-quasistatic behaviour 
of micro-refrigerators in underdamped as well as in overdamped cases in 
detail. Finally we conclude by focusing on the major differences between 
micro and macro machines manifested by the model system. 


2. The Model 


We consider a Brownian particle of mass m, position x and velocity v 
moving in a medium with friction coefficient 7 . The particle is trapped in 
a harmonic potential whose stiffness constant k(t) is periodically changing 
with time period r. The underdamped equation of motion of the particle in 
contact with a thermal bath of temperature T is given by mm 

mii = —qu — k{t)x + 7 / 7 T£(£), (1) 


where fluctuation dissipation relation between noise strength, temperature 
of the bath and friction coefficient is maintained. Throughout the paper 
Boltzmann constant ks is set to unity. In overdamped limit, the equation of 
motion reduces to 

7 x =-k(t)x + (2) 

The noises from the bath £ are Gaussian distributed with zero mean and 
are delta correlated, i.e., (£(£)) = 0 and (£(£i)£(t 2 )) = 25 (£1 — t 2 ). The time 
dependent stiffness of the trap is given by, 


m = a (l+l) 


= a/2 



= a /2 


0 < t < t/2 
t = t /2 
t/2 < t < T 

t = T. 


(3) 
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In the first step, the system undergoes an isothermal compression while in 
contact with the hot bath of temperature Th for 0 < t < r/2. During this 
step, work is done on the system. At t — r/2, the stiffness is instantaneously 
changed from a to a/2 (adiabatic expansion) and simultaneously the system 
is connected to the cold bath at temperature 7], disconnecting it from hot 
bath. In the third step, isothermal expansion is carried out where the stiffness 
is changed from a/2 to a/4 during r/2 < t < r. During this step work is 
extracted from the system. Finally the stiffness is changed by adiabatic 
compression from a/4 to a/2. In this step the system is detached from the 
cold bath and coupled with the hot bath again. A schematic diagram of the 
protocol is shown in Fig. [2]A. 


Isothermal process 
with hot bath at T/ L 



Adiabatic jump by 
sudden compression 
with switching from 
cold to hot bath. 


t 



Isothermal process 
with cold bath at T) 


Adiabatic jump by 
sudden expansion 
with switching from 
hot to cold bath. 


A 


B 


Figure 2: (Color online) A. Carnot-type refrigerating protocol: The red and grey lines de¬ 
pict the compression and expansion of the harmonic trap during the isothermal steps with 
temperature Th and 7] respectively. Black lines are the jumps of the protocol during the 
adiabatic steps. B. Schematic representation for a cyclic process of stochastic refrigerator 
operating between two reservoirs kept at temperatures Th and TJ. The cycle consists of 
two isothermal steps and two adiabatic steps according to the time varying protocol k(t). 
The violet line denotes a one dimensional potential V(x,t) and the filled region denotes 
the corresponding steady state distribution. 


One should note here that instantaneous jumps are not the only way 
to implement adiabaticity. There are examples of smooth protocols where 
one changes the temperature of the bath and the stiffness of the harmonic 
trap smoothly over time such that at long time limit the phase space volume 
remains constant, providing no average heat dissipation and consequently es¬ 
tablishing the adiabatic stepsjSTj. In this paper we restrict ourselves to this 
Carnot-type refrigerator protocol. This is the reverse of the protocol that 
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has been studied in |3Bj in context of single particle stochastic heat engine. 
There are other important protocols, for example Stirling-type protocol used 
in earlier studies |19|, where the adiabatic steps are absent. The compar¬ 
ative study between stochastic thermodynamics of micro machines driven 
by different protocols will be focused elsewhere. Our current work concerns 
stochastic thermodynamics of time periodic steady state (TPSS) of the sys¬ 
tem. It is well known fact [55] that when a Langevin system is subjected 
to a time periodic force, after initial transient (duration of which depends 
on the details of system parameter), the system will settle to a TPSS. In 
TPSS, the joint phase space distribution P ss (x,v,t ) is time periodic, i.e., 
P ss (x,v,t ) = P ss (x,v,t + r). 

3. Stochastic Energetics of the Single Particle Refrigerator 

Before investigating further, we briefly recapitulate the principle of stochas¬ 
tic thermodynamics. Using Eq. [lj change in internal energy, Au, for an 
underdamped Langevin system can be written as 



(4) 


w — q 


where u = \mv 2 + \k(t)x 2 , w — f Qjf cit and q = — f (—yu + y/yTQndt. 
The integrals in the second step are performed according to Stratonovich’s 
rule. Along a fluctuating trajectory of the particle, identifying the random 
variables u, w and q respectively as total internal energy of the particle, ther¬ 
modynamic work done on the particle and heat dissipated to the bath, one 
obtains first law of stochastic thermodynamics, valid for a single trajectory 
(Eq. g. Note that while w and q depend on the entire trajectory of the 
particle, Au depends only on its initial and final points. Similarly we can 
also write the first law in case of overdamped dynamics where u = \k{t)x 2 . 

For the system concerned here, work done on the particle in isothermal 
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steps for a typical cycle in TPSS is given by, 


^isoth 



dt 

dt 


du r , du 

dt + X /2 t dt 




( 5 ) 


The first integration in the r.li.s of above expression is defined along the path 
of isothermal compression and the second one is along isothermal expansion. 
By definition, in adiabatic steps q = 0 and therefore work in adiabatic steps 
during this cycle in TPSS is 


Wad 


A u 


r 

7 r_ M 


Ml - 

— u - 

+ 

L V2/ 

\2 J\ 




( 6 ) 


The first term of the above expression is the work done on the system along 
the path of adiabatic expansion while the second one is that along adia¬ 
batic compression. Simulating the dynamics of the particle via Eq. [lj we 
numerically calculate Wi SOt h and w a d to obtain the total work for a cycle as, 


W Wisoth T W a d- 


(7) 


Running the dynamics for a large number of cycles (N) in TPSS, we calculate 
work, averaged over all cycles, 

W = £ w ( 8 ) 

all cycles 


Using Eqj4]we calculate heat transfers in isothermal steps between the par¬ 
ticle and the baths ( q ). Heat transfer along isothermal compression is given 
by 


<h - - 


rr/2 . 

dt- ( kx 2 


T=T h 


+ 


u[ ^ ) - u[ 0 


and that along isothermal expansion is given by, 


(9) 





92 — — 


( 10 ) 










Since the heat transfers in adiabatic steps are zero, the total heat transfer in 
a cycle is 

q = qi + « 2 - (if) 

As before, running the cycle repeatedly for N times, we define the average 
heats, 

a = v 5Z qi] 02 = v 0 = v S q - (12) 

all cycles all cycles all cycles 

These non-quasistatic results should asymptotically match with quasistatic 
results which can be recovered as lirn^oo Q 3 = Q 3 with j = 1,2 and 
lim T _ >00 W = W t0 f Here QjS are the heat exchanges and W to t is the to¬ 
tal work in a cycle in quasistatic limit. Similarly, one can also calculate the 
change of total energy of the particle in a cycle, A u and its average over N 
cycles followed by its quasistatic counterpart U , taking r —* oo limit. There¬ 
fore, for the system concerned here, all the relevant thermodynamic variables 
involved in the first law of stochastic energetics and their path-averages, both 
in non-quasistatic as well as quasistatic regime, can be calculated by simulat¬ 
ing the underdamped and overdamped dynamics of the system. Quasistatic 
behavior of our system serves as an important benchmark for the simulation 
as it can also be obtained analytically. 

For any machine, quantification of its performance (e.g., efficiency) is one 
of the important issues. In case of microscopic engine [38], it has been shown 
that, in quasistatic regime average efficiency of the micro heat engine depends 
only on Similarly, in case of micro refrigerator we will determine average 
COP in quasistatic regime. We define COP for a cycle as 


w 

for a trajectory of the system. It fluctuates randomly from cycle to cycle. 
For large number of cycles, average stochastic COP is defined as, 

« = v £ e 

all cycles 

Due to fluctuations of q and w in both non-quasistatic and quasistatic regime, 
(e) 7 ^ = e, where e is the conventional definition of average COP. For 

finite time cycle, using fluctuation theorems [52, 53j IMj, it can be shown 


(14) 
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that e < e c = T ^/_ Ti ■ Note that this is valid for all cycle times. However, no 
such bound exists for (e). Being equipped with stochastic thermodynamics 
of our system, in the following sections, we explore quasistatic as well as 
non-quasistatic behavior of our system. 


4. Underdamped dynamics 

4.1. Quasistatic Limit 

Now we calculate thermodynamic quantities like, average work, heat ex¬ 
changes and internal energy changes for different sub steps of a cycle in 
quasistatic limit. In this limit, the duration of the protocol is longer than 
all other time scales. During isothermal processes, though the protocol is 
being changed, the system instantaneously adjusts itself to the equilibrium 
corresponding to the value of the protocol at that instant. Hence the work 
done along any isothermal process is the free energy difference between the 
final and initial state. In the first step, i.e., in isothermal compression, the 
work done on the system connected to the hot bath will be 


= n HP -> (7-/2)-) 
2 k(t = 0) 



(15) 


At t — ( t / 2 ) , the system is in equilibrium with hot bath and the corre¬ 
sponding distribution is given by 


Pr (x, v) = N\ exp 



(16) 


Here, N\ = is the normalization constant. The second step being instan¬ 
taneous, there will be no heat dissipation and the work done on the particle 
is the instantaneous change of its internal energy, given by 

W 2 = / \ (i ~~ °) x ' 2jP 2 v ) dxdv = -y- (!7) 


Similar to the first step, the work done on the particle in the third step is, 


W 3 


Ti k{t y r ) 
2 k(t = t/2) 



(18) 
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At t = t , the system will be in equilibrium with the cold bath. Hence, the 
corresponding distribution will be 


P T (x, v ) = N 2 exp 



(19) 


with N 2 = 4 ^?- Therefore, in the last step, i.e., in adiabatic compression 
process the average work done on the particle is 

W4 = / \ (i ~ l) x 2 p r( x ’ v ) dxdv = y • ( 20 ) 

The average work done during the full cycle in the quasistatic process is 


W tot =W 1 + W 2 + W 3 + W 4 

= ~T + ? + \ i ' Th ~ Ti)l ‘ l2 ' ( 21 ) 

As Th > Ti > 0, we find that W to t is always positive for any temperature 
difference. This implies that, in the quasistatic limit, on an average, work is 
always done on the system. 

To obtain heat absorption from the cold bath (Q 2 ) first we have to calcu¬ 
late internal energy change along the third step. The average internal energy 
at t = is U (| + ) = f Pz(x,v)dxdv = Since the system 

is in equilibrium with the cold bath at t = t~ , the average internal energy 
will be Ti. This leads to the change in internal energy in the third step, 
(Ti - ^). Using first law (Eq. 0 we obtain average heat dissipated to the 
cold bath, 

rri O rn 

Q 2 = — yin 2 — XJ -| (22) 

Similarly we can obtain the average heat, transferred to hot bath, 

T* ST 1 

Qi = -f ln2 ~ Th + ^ L . ( 23 ) 

Though Wtot is always positive, Q± and Q 2 can take negative as well as 
positive values depending on the ratio between hot and cold bath temper¬ 
atures. When Pr < 1.80, Q 2 is negative and is positive elsewhere. When 

< 2.29, Qi is positive and otherwise negative. The system will act as a 
refrigerator only when Q 2 < 0 and Q\ > 0 i.e., ^ < 1.80 in the quasistatic 
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limit. For 1.80 < jr < 2.29 system will act as a stochastic heater I (Fig. 
1) where both Q i and Q 2 are positive. When || > 2.29 heat flows from 
hot bath to the system (Qi < 0), but the system releases heat to the cold 
bath thereby acting as a stochastic heater of type 11. It is evident that even 
after using a Carnot-type refrigerating protocol one can obtain single particle 
refrigerator as well as heaters depending on Th and X) in quasistatic limit. 
Finally the COP of the system as a refrigerator in quasistatic limit is, 


-Q '2 T] In 2 + 27} 


3 T h 
2 


w t 


tot 


* _ L _ z _ 

f+T, + (T h -T,)ln2' 


(24) 


It should be noted here that even in the quasistatic limit, COP in our system 


is much smaller than the Carnot limit (e r = 


T h —Ti 


as shown in Fig. 



Figure 3: Comparison between quasistatic COP for our model in underdamped limit (e g , 
in red line) and the corresponding Carnot bound (e c , in green line) as a function of 
Blue filled circles denotes the values (obtained from numerical simulations) of e at different 
values of ^ for large cycle time. 


We now turn our attention to the results obtained by simulation in non¬ 
quasistatic regime. 

4-2. Numerical Results And Discussions 

We evolve the system using discretised Langevin dynamics with time 
step dt = 0.0002 in the underdamped as well as overdamped limit [Eq. [I] 
and Eq. [ 2 ]. The system is driven by time periodic protocol [Eq. [ 3 ]. We 
follow Heun’s method [58]. We have set 7 = 1 and rri — 1. All the physical 
quantities are in dimensionless form. Throughout the paper we have fixed 
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a = 5 and T) = 0.1. We make sure that, after the initial transient regime 
10 3 cycle time), the system settles to a TPSS i.e., P ss (x, v, t+r ) = P ss {x, v, t). 
We use the same discretisation to numerically calculate the thermodynamic 
variables along a trajectory. We consider at least 10 5 cycles of operations 
and thermodynamical quantities are averaged over all these cycles. 

Phase diagram: For each (r, T h ) pair, we calculate W, Q i, and Q 2 with 
N ~ 10 5 and a phase diagram [Fig. [4] for the operational modes of the 
system is obtained. From the phase diagram it is evident that when r < 2, 
for any Th, Qi < 0, Q 2 > 0 and W > 0. In this mode, on an average the 
system heats up the cold bath. Therefore it is a stochastic heater of type II. 


1 


T h 


o.i 

071 1 10 100 

r 

Figure 4: Phase diagram for underdamped motion: Different modes of operation of our 
system following underdamped Langevin dynamics with Carnot refrigerating protocol. 
Open boxes(pink): refrigerator, asterisk(blue): heater-I, cross(green): heater-II, plus(red): 
engine. Here all the points in the phase diagram are obtained by averaging over 10 5 
trajectories. 



For large r(~ 100) and high T h , Qi and W goes into the system whereas 
Q 2 comes out of the system to heat up the cold bath, providing us again 
heater-II. At large r, but with intermediate Th, we obtain heater-I where 
Qi flips its direction of flow to heat up the hot bath keeping everything else 
the same as before. If we reduce Th further, we obtain stochastic refrigerator 
where, even Q 2 flips its direction and flows from cold bath to the system while 
others remain same as heater-I. Interestingly for Th > 0.35 and 2 < r < 10, 
we observe that on an average, the system can do work, taking heat from 
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the hot bath and releasing it partially to the cold bath. It implies that 
even if we use a Carnot-type refrigerating protocol, deep in non-quasistatic 
regime, for a particular range of r and Ty ll our system can behave as a 
stochastic heat engine. Therefore, if we consider both non-quasistatic and 
quasistatic regimes, the system can act in all thermodynamically possible 
modes, depending on the values of (r, T h ). It should be mentioned here that 
for large r, phase boundaries are correctly predicted by the quasistatic limits 
calculated analytically. The details of the phases and phase boundaries are 
depicted in Fig. [4j Note that, in the phase diagram crossing any phase 
boundary once, implies flipping the direction of any one of the Q\. Q 2 and 
W. 

One can further explore different modes of operations in Fig. [5] where we 
have plotted W, Q\ and Q 2 with respect to cycle time r, for different Th. If 
temperature difference between two reservoirs is small, W is always positive 
and Qi is negative for small r but positive for large r. The behavior of Q 2 
with t is opposite to that of Qi. Thus for small temperature differences, if we 
vary r from small to large values, we obtain heater-II first and then heater-I 
and finally a refrigerator. For higher temperature differences, with increasing 
r, Q\ becomes negative and Q 2 becomes positive. W shows a remarkable 
non-monotonic (convex) behavior within a range of r where it shows a dip 
that can even be below zero for higher temperature differences. Therefore, 
in this regime of temperatures, though we obtain heaters with very small 
and large r, our system acts as a stochastic heat engine within a particular 
range of intermediate r. W, Q\ and Q 2 ultimately saturates with r to values 
as found in quasistatic calculation. We note from Fig. [5]4 that W becomes 
independent of Th for a given X) at two different r within which W is non 
monotonic. 
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Figure 5: A: W versus r for different T hl B: Q 1 versus t, C: Q 2 versus r, D: A S tot versus 
r. Red: Th = 0.12, Green: Th = 0.2: , Blue: Th = 0.5, Cyan: Th = 1.0). Two crossing 
points of all the curves in A are marked by broken black circles in A. At these points W 
is independent of Th¬ 


ill Fig. HF we have plotted average total entropy production (A S to t = 
j r + jr) per period. It shows that initially AS tot increases with r. Above 
r = 10, A Stot tends to saturate to quasistatic value. This sets a characteristic 
time scale in our system. The characteristic time scale that determine the 
relaxation time of our system is r ~ 10. The entropy production increases 
monotonically with T h . In Fig. [5]A, we see two special points where the value 
of W for particular values of r are same even though the hot bath temperature 
T h is changed. We explore the probability distribution P(w) of work at those 
two values of r for different Th as shown in Fig. [GjV and [6)3. It is clear from 
the distributions that the fluctuations increases with increasing Th, keeping 
the mean unchanged. P(w) is asymmetric with positive mean. Finite weight 
for w < 0 arises from realisations which doesn’t act as refrigerator. 
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Figure 6: Work distributions are plotted for various bath temperatures : T h = 0.12 (red), 
Th = 0.2 (green), T), = 0.5 (blue), Th = 1.0 (cyan). The cycle time for A and B are 
r = 1.916 (first crossing point) and r = 14.98 (second crossing point) respectively. The 
distributions are obtained by considering 10 5 trajectories. 


Refrigerator mode: In the refrigerator regime of the phase diagram, most 
of the trajectories followed by the system maintains the refrigerating condi¬ 
tion: qi > 0, q 2 < 0 and w > 0. However, there are considerable number of 
trajectories where the system does not follow this condition. This number 
reduces with increasing r but still remains finite even in quasistatic regime. 
Due to fluctuations, the difference between (e) and e (Fig. [7] A) is prominent 
even at large r. The saturation values of (e) and e (as given in TABLE I) 
are far below than the Carnot limit (e c = 5 in our case). It should be noted 
that the saturation value of e matches with our quasistatic result. 



Figure 7: Variation of (e) (black, broken) and e (red, continuous) with r at 7^=0.12. The 
average quantities are obtained by considering 10 5 trajectories. 


P(e), the distribution of stochastic COP, is shown in Fig. [8] for two 


16 














different cycle times, r = 10 and r = 100 (both for T\ = 0.12). e can take 
values from — oo to oo. There are considerable number of realizations where 
e can be negative or can occur beyond e c . Moreover, we have noticed (see 
inset of Fig. [8]) that the tail of P(e) decays as a power law ( e“) for several 
decades. The exponent a depends on the system and protocol parameters 
(e.g., r) and are given in figure captions. So far to our knowledge, no study 
exist for the probability distribution and its power law tails for stochastic 
COP. It is not clear whether this power extends indefinitely for large values 
of e. However, given our numerical data we can calculate variance which is 
finite. 

For both distributions, variance a e > (e) (see TABLE [I] for more details). 
Thus (e) ceases to be a good physical variable and one has to resort to the 
study of full probability distribution. From Fig. [8] it is also apparent that, 
P(e) becomes sharper with increasing r. From the phase diagram we know 
that, with r ~ 10 2 , our system is closer to the quasistatic limit. Therefore, 
even in quasistatic regime, fluctuations in e are significant. 



Figure 8: Probability distribution of e at r = 10.0 (red, continuous) and r = 100.0 (black, 
broken) at T h = 0.12. Inset shows positive tails of the distributions plotted in log scale. 
They behave as e“. For r = 10.0 a =-1.897 ± 0.004 ( slope of the red, continuous line) 
and for r = 100.0 a = -1.875 ± 0.004( slope of the black, broken line). 


In Fig. [9] distribution of COP is plotted at r = 100 for different numbers 
of cycle(n). As n increases, distribution becomes sharper. The tail on the 
negative side disappears thereby indicating the fact that realisations leading 
to non-refrigerator mode reduces drastically. Hence our system works as 
a reliable refrigerator in large n limit. Tail on the positive side also gets 
suppressed as n increases. The standard deviation decreases as n increases. 
For n > 50, standard deviation becomes less than the mean value. This 
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implies that for large n, (e) becomes a well defined physical quantity. Average 
COP does not approach Carnot limit for large n. 



Figure 9: Probability distribution of e for different cycle numbers(n) at r = 100. 


In Fig. [To] A and Fig. [Tb) f we have plotted P(w, 1 / 2 ), the joint probability 
distributions of w and q 2 over 10 5 trajectories, at T h = 0.12 for two different 
cycle times r = 10 and 100. For a particular cycle, e can be positive when q 2 
and w has opposite signs, but only fourth quadrant of these joint probability 
distribution represent a refrigerator operation. In Fig. 0c only 46.2% tra¬ 
jectories follow the refrigerating condition whereas in Fig. [Tofrh it is enhanced 
to 52.1%. Therefore, the system will behave more reliably as a refrigerator 
if we run it close to quasistatic limit. Table I summarizes important results. 



Figure 10: Joint distributions: A: Distribution of w and q 2 at r = 10.0, Th = 0.12, B: 
Distribution of w and q 2 at r = 100.0, T h = 0.12. 


Engine mode: From the phase diagram, it is apparent that the system 
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Table 1: 


T 

w 

Qi 

Q\ 

e 

<*> 


acts as a 
refrigerator 

10.0 

0.0277 

0.0527 

-0.0250 

0.903 

0.458 

10.3 

46.2 % 

100.0 

0.0271 

0.0696 

-0.0423 

1.561 

0.929 

11.7 

52.1 % 


acts as an engine only in the non-quasistatic regime where fluctuations are 
dominant. The joint distribution P(w,qi) at (r = 5.0 ,T h = 1.0), a typical 
point within engine region of the phase diagram, shows that for only 47.3% 
of ~ 10 5 trajectories, the condition for engine, i.e. q\ < 0, q 2 > 0 and w < 0, 
is followed (Fig. 0 - 



w 

Figure 11: Joint distribution of w and q± at r = 5.0,77 = 1.0. 


The effect of fluctuation is also reflected in the distribution of stochastic 
efficiency r] (Fig. 12). Thus this engine is unreliable. We notice that the 


efficiency distribution is bi-modal. Both the peaks as well as (r/) are well 
below the Carnot bound q c = 0.9. Due to strong effect of fluctuations, 
the standard deviation of efficiency a v , is much larger than its average (?/). 
The tail of P(r/) decays as a power law ( rj 13 ) for several decades with (3 = 
—2.143 ± 0.010. The exponent is different from that we have obtained in 
case of P(e). It suggests that altering r and Th not only implies the change 
of the thermodynamic mode of operation in r — Th space but it also implies 
significant difference in the behaviour of probability distributions of their 
performance index (efficiency / COP). 
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1.2 



Figure 12: Probability distribution of rj at r = 5.0 at T ^ = 1. Here ( 77 ) = 0.24 and 
a v = 3.83. Inset: positive tail of the distribution, plotted in logarithmic scale, shows a 
power law decay with exponent -2.143 ± 0.010 (slope of the straight line). 


5. Overdamped dynamics 

The difference between the overdamped and nnderdamped dynamics is 
the absence of inertia. The relevant equation of motion for overdamped 
dynamics of the particle Eq. [2] has no acceleration term. This approximation 
is valid when the time steps of the observation are much larger than mj 7 . We 
now discuss the performance of our system following overdamped dynamics. 


5.1. Quasistatic Limit 

In the quasistatic limit the calculation of thermodynamic quantities are 
similar to that of the underdamped case. The main difference is that, in 
absence of inertia, only potential energy contributes to the internal energy. 


The work done along adiabatic steps will be same as given by Eq. [17] and 
Eq. [20| except for the fact that the probability distribution will now depend 
only on position of the particle not on the velocity. The work done in the 
isothermal process will be exactly same as in Eq. [15] and Eq. 18 The 


expression for total work done on the system in this process will be the same 
as that we obtained in the underdamped case (Eq. |2ll) 


At t = 


T + 
2 ’ 


average internal energy is t/(| + ) = -p. At t — t system is 


in equilibrium with the cold bath and corresponding average internal energy 
is U(t~ ) = y. The average internal energy change in isothermal expansion 
process is (y — y). Hence, the average heat that goes to the cold bath in 
quasistatic limit is given by 


Qo = — — In 2 — — + —. 
^ 2 2 4 


( 25 ) 
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The expression for Q i, the average heat that is exchanged between hot bath 
and the system, is 

Qi = yln2 — y + 7], (26) 

For overdamped case, the system will act as a refrigerator when < 3.39. 
For 3.39 < ^ < 6.51, the system will act as stochastic heater-I. Beyond this 
limit we get stochastic heater-II. The coefficient of performance for stochastic 
refrigerator is given by 


69 Wi 


-Q 2 Ti In 2 + T; —y 

tot — + T) + (T/,, — Tj) In 2 


(27) 


which is again below the Carnot bound as shown in Fig. 13 



Figure 13: Comparison between quasistatic COP for our model in overdamped limit (e g , 
in red line) and the corresponding Carnot bound (e c , in green line) as a function of . 
At T rp = 2, e q equals Carnot limit. Blue filled circles denotes the values (obtained from 
numerical simulations) of e at different values of for large cycle time. 


We would like to emphasize that only for jr = 2, e q equals Carnot limit. 
This is an interesting observation [5,9]. It is important to note that the 
adiabatic jumps, being instantaneous, the probability distribution remains 
unchanged during this process. The system has to relax to new equilibrium 
along the isotherms after sudden change in the temperature. The relaxation 
time is negligible in the quasistatic limit i.e. when the cycle time r is in¬ 
finitely large in comparison to the relaxation time. This relaxation leads to 
an additional heat flow which accounts for the change in the internal energy 
during the relaxation process. This additional heat flow becomes zero if the 
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system is in equilibrium immediately after the jump. This can be achieved 
either by changing the jump values of the protocol or by tuning the tem¬ 
peratures. In the overdamped case, we have found that for ^ = 2, this 
additional heat flow is zero and hence the Carnot bound is obtained for our 
specified protocol. However, in the underdamped limit, one cannot make 
this additional heat flow equal to zero for any value of adiabatic jump or jr 
and we do not acheive Carnot bound in this quasistatic limit. We have also 
found for engine protocol [38] that Carnot efficiency e c = 1 — ^ is achieved 
only for the value ^ = 2 in the overdamped limit. Thus, in the overdamped 
limit, refrigerators/engines work in reversible mode only for values ^ = 2 
for our protocol. We will show later that for this particular value, total 
entropy production over a period vanishes. We now consider the results in 
non-quasistatic regime obtained in our simulations. 


5.2. Numerical Results and Discussions 

We have scanned the parameter space ( r — T h ), keeping Ti fixed at 0.1, to 
obtain the phase diagram in overdamped regime (Fig. 14). It clearly depicts 
four different modes of operation of the system, as described in Fig. [l} In 
contrast to the underdamped case, no critical cycle time is required for the 
operation in refrigerator mode. The system acts as a refrigerator for higher 
temperature differences compared to the earlier case. Therefore, total phase 
space area of the refrigerator mode has increased in this limit. The phase 
boundaries in quasistatic limit are consistent with our analytical results. The 
engine region is clearly visible for 0.02 < r < 1.0 and Th > 0.3. 
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Figure 14: Phase diagram for overdamped motion: Different modes of operation of our sys¬ 
tem following overdamped Langevin dynamics with Carnot refrigerating protocol. Open 
boxes(pink): refrigerator, asterisk(blue): heater-I, cross(green): heater-II, plus(red): en¬ 
gine. The averaging, considered here, is over 10 5 trajectories. 


To explore further, in Fig. [15] we have plotted W, Q\ and Q 2 as a 
function of r for different T Similar to the inertial case, W is positive for 
all r at low T h (— 0.12,0.2). ffowever at higher values of ^(=0.5, 1.0), non¬ 
monotonicity of W with respect to r becomes apparent. For smaller r it 
is negative, showing a dip around r ~ 0.8. After the dip, W increases and 
eventually saturates for large r at a positive value, predicted by quasistatic 
result. Interestingly, in high friction limit, the range of r where W is non¬ 
monotonic, starts from r ~ 0 whereas in the underdamped case, the range 
was bounded between two nonzero cycle times. From Fig. 15 it is apparent 
that behavior of W, Q\ and Q 2 with r, are all consistent with the phase 
diagram for different modes of our system in the high friction limit. 

In Fig. 15 3 we have plotted A S tot = f 1 + fr as a function of r. The 


characteristic time scale that determine the relaxation time of our system in 
overdamped limit is r ~ 5. For r > 5, A S tot tends to saturate to quasistatic 
value. Interestingly, for ^ = 2, A S to t approaches to zero in the quasistatic 
limit. It implies, for this special case, refrigerator works in a reversible mode 
and the corresponding COP equals Carnot value as discussed earlier. For 
all other values of A., our system works in an irreversible mode with finite 
value of A S to t even in quasistatic limit therby preventing it to reach Carnot 
value. In Fig. [15}\, we see that at r = 1.957, W does not change with Th- 
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Figure 15: In high friction limit, A: W versus r for different Th, B: Qi versus r, C: Q 2 
versus r, D: (A S to t) vs r. Red: = 0.12, Green: Th = 0.2: , Blue: Th = 0.5, Cyan: 

Th = 1.0). W becomes independent of Th at the crossing points of all curves in A, as 
shown within broken black circle. Here we take the average over 10 5 trajectories. 


In Fig. 16 we have plotted P(w) for different to show that the mean is 
constant but the width increases with Th . P(w) is asymmetric with positive 
mean. Finite weight for w < 0 comes from realisations which doesn’t act as 
refrigerator. 



Figure 16: Work distributions are plotted at r = 1.957 i.e., at the crossing point of the 
W vs. r plot, for various bath temperatures : Th = 0.12 (red), Th = 0.2 (green), Th = 0.5 
(blue), T h = 1.0 (cyan). 
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Figure 17: Probability distribution of e at r = 10.0 (red, continuous), t = 50.0 (blue, 
broken) and r = 100.0 (green, broken) at Th = 0.12. Inset:positive tails of the distributions 
plotted in log scale. Slopes of the straight lines indicate the exponents. For r = 10.0, 50.0 
and 100.0, a =-2.255 ± 0.010 (red, continuous line), -2.275 ± 0.010, -2.288 ± 0.008 (green, 
broken line) respectively. 


Iii Fig. 17 we have plotted the distribution of e, P(e), for r = 10, 50 and 
100. Close to quasistatic limit (r = 100) we observe a dip in P(e) at e = 1 and 
the distribution becomes bi-modal. The fluctuations prevail even at large r, 
thereby making e a non-self-averaging quantity. For large e, P(e) shows power 
law tails with exponents given in figure caption. Notably exponents are lesser 
than —2, indicating finite variance. The exponent values are different from 
that obtained in the underdamped case. This suggests that the exponents 
are not universal as our system works in an irreversible mode. Mean and 
standard deviations of different stochastic thermodynamical variables for two 
different values of r at Th = 0.12 are given in Table 2. 


Table 2: 


T 

w 

Qi 

Q'2 

e 

<e) 

cr e 

acts as a 
Refrigerator 

10.0 

0.0272 

0.0766 

-0.0494 

1.830 

0.916 

6.863 

65.6 % 

100.0 

0.0270 

0.0812 

-0.0541 

2.004 

1.348 

7.80 

67.1 % 


We have plotted the distribution of stochastic efficiency in Fig. [18] at 
r = 0.5 and Th = 1.0, a typical point of the phase diagram where our 
system works as stochastic heat engine. We found that only for 40.3% cycles, 
the system works as a heat engine. Thus, even in overdamped regime, the 
engine operation is unreliable. Both P(e) and P(r/) show power-law decays 
in their tails for several decades. The exponent values are given in the figure 
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captions. The exponent of P(e) is different from that we have obtained in 
case of underdamped dynamics. 



Figure 18: Probability distribution of rj at r = 0.5 at T h = 1. Here ( 77 ) = 0.54 and 
ay = 3.56. Inset: the positive tail of the distribution behaves as a power law with exponent 
-2.155 ± 0.014 ( slope of the straight line). 


6. Conclusion 

We will conclude by focusing on the main results to point out the crucial 
differences between macro and micro refrigerators. In the macro-world, ide¬ 
alized refrigerator is like an engine running backward in time, i.e., using the 
work to extract heat from the cold bath and dumping the heat at hot bath. 
We have shown that such a simple picture is not valid for single particle heat 
engine at nano-scale. The engine running in a reverse order may produce 
refrigerator, stochastic heater of type I and II as well as unreliable heat en¬ 
gine. For our model of study, fluctuations in stochastic COP and efficiency 
dominate their mean values even in the quasistatic case. This implies that 
in such a situation mean is not a good physical variable and one must study 
the behaviour of full probability distributions which in all our cases contain 
power law behaviour in their tails with varied exponents. Stochastic COP 
and efficiency may exceed Carnot bound. However, averaged COP and effi¬ 
ciency, as defined here, are bounded by the Carnot value. This can readily 
be shown by fluctuation theorem of heat engines. 

In every cycle, a macro thermal machine under refrigerating protocol acts 
as a refrigerator. This is not true for micro refrigerators. We know that the 
trajectory dependent work, heat and COP are randomly distributed variables 
over the cycles of a micro refrigerator. The cycle time can be large or small, 
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but there will always be a considerable number of cycles where the system 
will not run as a refrigerator due to fluctuations of these thermodynamic 
variables. This makes the micro refrigerator unreliable in comparison to 
macro refrigerators. This reliability, quantified as a fraction of cycles running 
as a refrigerator, increases as we increase the cycle time r. The obtained heat 
engine is most unreliable. 

The phase diagram and thermodynamic quantities describing the micro 
refrigerators as well as the micro heat engines are crucially protocol depen¬ 
dent. We have noticed in our other studies (will be published elsewhere), 
that the micro refrigerators / heat engines, running under Stirling protocol 
(which is very similar to Carnot protocol but devoid of the adiabatic jumps), 
though produce qualitatively similar phase diagrams as Carnot-protocol, are 
less fluctuating in large cycle time limit. Moreover, only in the underdamped 
case in a small parameter space heat engine operation is possible. The adi¬ 
abatic jumps are key players behind the unreliability of micro heat engines 
in [3E] or refrigerators described in earlier sections. Only in the overdamped 
quasistatic limit we have shown that for our given specific protocol, refriger¬ 
ator can work in a reversible mode of operation provided = 2. This clear 
from zero entropy production for this case (e.g., see Fig. [l5]D). Power output 
is found to be an useful optimization criterion for finite time heat engines. On 
the other hand, many optimization criteria were proposed in case of refrig¬ 
erator. Widely used target function was introduced in [HU ET, S2) ESI iMj by 
giving equal footing for COP and cooling rate, known as ^-criterion [ 621, [63] . 
In view of new findings in Brownian refrigerators, it will be of immense inter¬ 
est to study these quantities in our system. It may be noted that our studied 
engine or refrigerator do not reach Carnot limit in the efficiency / COP in 
the quasistatic regime. The engine is not a microscopic version of Carnot 
engine. There are other engines which reproduces Carnot limit for efficiency 
in the quasistatic regime [20]. Studies on differences in these two types of 
engines will be useful. Work along this direction is in progress. 
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